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ABSTRACT 

Context. Circumstellar disks and outflows play a fundamental role in star formation. Infrared spectro-interferometry 
allows the inner accretion-ejection region to be resolved. 

Aims. We study the disk and Br7-emitting region of MWC 297 with high spatial and spectral resolution and compare 
our observations with disk-wind models. 

Methods. We measured interferometric visibilities, wavelength-differential phases, and closure phases of MWC 297 with 
a spectral resolution of 12 000. To interpret our MWC 297 observations, we employed disk-wind models. 
Results. The measured continuum visibilities confirm previous results that the continuum-emitting region of MWC 297 is 
remarkably compact. We derive a continuum ring-fit radius of ~2.2 mas (~0.56 AU at a distance of 250 pc), which is ~5.4 
times smaller than the 3 AU dust sublimation radius expected for silicate grains (in the absence of radiation-shielding 
material). The strongly wavelength-dependent and asymmetric Br7-emitting region is more extended (~2. 7 times) than 
the continuum-emitting region. At the center of the Br7 line, we derive a Gaussian fit radius of ~6.3 mas HWHM 
(~1.6 AU). To interpret the observations, we employ a magneto-centrifugally driven disk-wind model consisting of an 
accretion disk, which emits the observed continuum radiation, and a disk wind, which emits the Br7 line. The calculated 
wavelength-dependent model intensity distributions and Br7 line profiles are compared with the observations (i.e., 
A"-band spectrum, visibilities, differential phases, and closure phases). The closest fitting model predicts a continuum- 
emitting disk with an inner radius of ~0.3 AU and a disk wind ejection region with an inner radius of ~0.5 AU 
(~17.5 stellar radii). We obtain a disk- wind half-opening angle (the angle between the rotation axis and the innermost 
streamline of the disk wind) of ~80°, which is larger than in T Tau models, and a disk inclination angle of ~20° (i.e., 
almost pole-on). 

Conclusions. Our observations with a spectral resolution of 12 000 allow us to study the AU-scale environment of 
MWC 297 in ~10 different spectral channels across the Br7 emission line. We show that the Jf-band flux, visibilities, 
and remarkably strong phases can be explained by the employed magneto-centrifugally driven disk wind model. 

Key words. Stars: individual: MWC 297, Stars: pre-main sequence, Stars: winds, outflows, Stars: circumstellar matter, 
Techniques: interferometric, Techniques: spectroscopic 



1. Introduction 

Circumstellar disks and outflows play a fundamental role in 
star formation. Testing theoretical accretion-ejection mod- 
els is very challenging since important processes occur on 
AU and sub-AU scales, which cannot be resolved by conven- 
tional imaging. However, infrared spectro-interferometric 
observations with milli-arcsecond spatial resolution and 
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* Based on observations made with ESO telescopes at Paranal 
Observatory under programme ID 081.D-0230(A). 



high spectral resolution allow unprecedented studies of the 
inner accretion region, which hosts fascinating astrophysi- 
cal objects such as inner gaseous accretion disks, dust disks, 
disk winds, jets, etc. Interferometry in the near- and mid- 
infrared is able to probe the gas and dust distribution on 
these scales and constrain their physical properties. These 
observations allow us to study the nature of the accretion 
process and the launching of jets and winds. 

With a spectral typ e of B1.5V, a mass of ~10 Mq, and a 
distance of 250 ±50 pc (jDrew et alJll997ft . MWC 297 is one 
of the nearest massive young stars. The geometry of the disk 
of MWC 297 has been investigated by various techniques 
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at near-infrared (iMillan-Gabet et al. _ 2001 ; lEisner et al 



2004; 
2005 



Monnier et al.ll2006HMalbet et al.ll2007t lAcke et al 



2008) , and radio wavelengths ([Mannings fc Sargent 19971 : 



Dre w et al.lll997l ). Interestingly, several interferometric ob- 



serva t ions of MWC 297 (lEisner et al.ll2004bHMonnier et~al 



Acke et al.1 120081: iKraus et al 



2005; IMalbet et al.l 12007 

2008al) have shown that the compact characteristic size 
of the near-infrared continuum-emitting region is several 
times smaller than the theoretically expected ^3 AU 
radius of the dust sublimation rim (the exp ected radius if 
there is no shielding of the stellar radiation; [Monnier et al.l 
120051) . 

The estimates of the inclination angle of M WC 297's cir- 
cumst ellar disk have been highly controversial. iDrew et al.l 
(|l997l) argued for an almost edge-on disk orientation. 
Moderate inclination angles of the circumstellar disk of 
MWC 297 were found wi th infrared in terferometric ob- 
For instance, IMalbet et "all 



servations. hoi instance, IMalbet et al.l (|2007l) obtained 
VLTI/ AMBER observations with a medium spectral res- 
olution of 1500 and measured the wavel ength depen- 
dence of the visibility across the Br7 line. IMalbet et al.l 
(|2007h interpreted their spectro-interferometric measure- 
ments of MWC 297 with an anisotropic s tellar- wind model 
(jStee &: de Arauio 1994; St ee et aT1ll995l) and found an in- 
clination of ~20°. lAcke et all (|2008l ) obtained both near- 
infrared AMBER H- and if-band observations and mid- 
infrared MIDI observations and concluded that the inclina- 
tion is smaller than 40°. 

In this paper, we present the first spectro- 
interferometric measurements of the inner circumstellar 
environment of MWC 297 with AMBER's high spectral 
resolution mode (i?= 12 000). This high spectral resolution 
allows us to spectrally resolve the B1-7 emission line 
and measure the interferometric observables (spectrum, 
visibilities, wavelength-differential phases, and closure 
phases) in ^10 different spectral channels across the line. 
In Sect. 2, we present the AMBER observations. In Sect. 
3, we derive characteristic continuum and Br7 line sizes 
of the circumstellar environment using simple geometric 
models. In Sect. 4, we present radiative transfer modeling 
of the disk and disk wind and compare the if-band flux, 
visibilities, and phases with the model predictions in the 
continuum and several spectral channels across the Br7 
emission line. In the Appendix, we discuss the disk-wind 
modeling, spectroscopy of MWC 297, wavelength cali- 
bration, and the use of AMBER's Beam Commutation 
Device. 



2. Interferometric VLTI/AMBER/FINITO 

observations with spectral resolution of 12 000 

We observed MWC 297 on 6 April 2008 with ESO's Very 
Large Telescope Interferometer (VLTI) and its AMBER 
beam combiner instrument in the course of t he AMBER- 
GTO program 081.D-0230(A). AMBER (|Petrov et al.l 
l2007l) allows interferometric three-telescope observations 
with low, medium, or high spectral resolution. The obser- 
vations are described in Tabled] For these observations, the 
E0-G0-H0 array of the 1.8 m auxiliary telescopes (ATs) and 
AMBER's high resolution mode (HR mode; R = 12 000) 
were employed. To obtain interferograms with a high signal- 
to- noise ratio (SNR), we used the fringe tracker FINITO 
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Fig. 1. AMBER observations of MWC 297 with spectral 
resolution of 12 000: (from top to bottom) wavelength de- 
pendence of flux, visibilities, wavelength-differential phases 
(for better visibility, the differential phases of the shortest 
and longest baselines are shifted by +40 and —40°, respec- 
tively), and closure phase observed at projected baselines 
of 14.0, 28.0, and 42.1 m along the position angle (PA) of 
68.0° on the sky. The wavelength scale at the bottom is the 
directly observed one, i.e. without heliocentric or LSR cor- 
rection (see text). However, the radial velocity scale at the 
top, gives spectrum, visibilities, and phases as a function of 
the Br7 Doppler shift in the Local Standard of Rest (LSR) 
frame. The dashed vertical line indicates the centroid Br7 
vacuum wavelength (2.1661 \xm) in the LSR frame. 
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and a detector integration time of 8.0 seconds per inter- 
ferogram. With excellent atmospheric conditions (seeing 
0.4-0.6"), the FINITO instrument in the co-phasing mode 
allowed us to stabilize the atmospheric and mechanically 
induced fringe motion. FINITO records temporally modu- 
lated £f-band interferograms with a high frequency, from 
which the atmospheric phase shift is determined and then 
compensated for using the VLTI delay lines. 

For data reduction, we used AMBER's amdlib data 
reduction softwa re (release 2-2), which employs the 
P2VM algorithm (|Tatulli et alJl2007t ) to derive wavelength- 
dependent visibilities, wavelength-differential phases, and 
closure phases. In addition, the spectrum can be extracted 
as an additional observable. 

Along with the science observation of MWC 297, we ob- 
served the interferometric calibrator star HD 175 583, which 
allowed us to calibrate the atmospheric transfer function. 
We recorded 100 target interferograms (plus 40 calibrator 
interferograms) without AMBER's beam commutation de- 
vice (BCD OUT data) and 100 target interferograms plus 
40 calibrator interferograms with the BCD (BCD IN data). 
The BCD allows the study of systematic instrumental ef- 
fects (see Appendix [C]). A fraction of the interferograms 
were of lower quality because of the lower FINITO perfor- 
mance. Therefore, we selected the 50% of the BCD OUT 
data and 50% of the BCD IN data with the highest fringe 
SNR among both the target and the calibr ator interfer- 
ogram s to improve the visibility calibration (jTatulli et al.l 
2007). From this data, we derived the spectrum, visibili- 
ties, wavelength-differential phases, and closure phases, ob- 
served along the position angle (PA) of 68.0° on the sky (see 
Fig- Hi- The differential phases and closure phases derived 
from both independent data sets (BCD OUT and BCD IN) 
are shown in Fig. IC.ll of Appendix [Cj Subtraction of the 
closure phases obtained with and without BCD cancels po- 
tential instrumental phase drifts. Figure IC.2I presents the 
derived closure phase and a fit of a constant to the con- 
tinuum closure phases. The obtained averaged continuum 
(2.147-2.194 fj,m) closure phase is -1.28 ± 0.42°. 

With the high spectral resolution of R = 12 000, it is 
possible to spectrally resolve the Br7 line and measure 
the visibilities and phases in ~10 different spectral chan- 
nels across the Doppler-broadened Br7 line. The visibilities 
show a strong wavelength dependence of the size of the line- 
emitting region. The wavelength-differential phases and the 
closure phase indicate remarkably strong photocenter shifts 
and asymmetries, respectively, of the circumstellar environ- 
ment in many spectral channels across the Br7 line. We 
show in Sect. U] that this strong wavelength dependence and 
the large values of the phases can be explained by disk- wind 
models. 

The wavelength calibration of the AMBER data in Fig. 
[I] was accomplished using the many telluric lines present in 
the region 2.15-2.19 jum, as described in the Appendix iBl 
We estimate an uncertainty in the wavelength calibration 
of ~3 kms . The spectrum in Fig. [1] (top) was normal- 
ized to the nearby continuum. The wavelength scale at the 
bottom is shown without heliocentric or Local Standard of 
Rest (LSR) correction. However, the radial velocity scale at 
the top of Fig. [T]is calibrated to give spectrum, visibilities, 
and phases as a function of the Br7 Doppler shift in the 
Local Standard of Rest (LSR) frame (16.5 kms -1 has to 
be subtracted to convert the LSR values to the heliocentric 
frame). 



The obtained Br7 line profile has an equivalent width 
of -17.2 A. Using AMBER's HR mode, we were able to re- 
solve the Br7 line and determined a FWHM of ^60 kms -1 . 
By performing a Gaussian fitting of the line profile, a cen- 
troid LSR velocity of 13.8 ± 3 kms -1 was obtained. No 
P-Cygni absorption or double-peaked emission is seen. The 
FWHM of the Br7 line is compatible with the value of 
53 kms" 1 found by iDrew et all (|1997ft for Bra. On the 
other hand, the centroid velocity of Br7 is more similar to 
the one m easured for Ho? a nd H/3 (16 and 19 kms" 1 , re- 
spectively; IDrew et al.lll997t ) than to that of Bra (3 kms" 1 , 
IDrew et al.lll997t) . 



3. Comparison of the observations with geometric 
models 

3.1. Characteristic size of the continuum-emitting region 

To measure the characteristic size of the circumstellar en- 
vironment in the K-band continuum, we fitted geometric 
models to the visibilities obtained along the PA of 68.0° o n 
the sky (see Table©. As pointed out bv lAcke et all (|2008h . 
who observed MWC 297 with the same VLTI E0-G0-H0 ar- 
ray, one-component uniform-disk, ring, and Gaussian disk 
models are unable to fit the shape of the continuum visi- 
bility. Therefore, we fitted a two-component Gaussian plus 
a point sourc e with a 10% flux contribution (derived from 
the SED; as in lAcke et al.ll2008h to the wavelength-averaged 
continuum visibilities. Our observations are averages over 
the continuum wavelengths between 2.147 and 2.194 /im 
(except the line region) because of the narrower wavelength 
range of AMBER's HR mode. 

Figure [5] (left) shows that the fitted two-component 
Gaussian consists of a dominant compact Gaussian with 
a FWHM diameter of 4.63 ± 0.21 mas (corresponding to 
1.16 ±0.05 AU at a distance of 250 pc) in addition to an ex- 
tended Gaussian halo with a FWHM diameter of ~22 mas 
(see Table [5]). Owing to the lack of data at low spatial 
frequencies, the size and shape of the extended halo com- 
ponent is not well constrained. The flux ratio of the domi- 
nant compact component to the extended halo component 
is 4.6. Our FWHM diameter of ~4.6 mas obtained for the 
dominant compact Gaussian component is in good agree- 
ment w ith the diameter of ~4.3 mas reported bv lAcke et al.l 
d2008h . 

To verify MWC 297's location i n the size-luminosity 
relation (where rin g fits are used; iMonnier et al.l [2005; 
iKraus et alJl2008ah . we determined its ring-fit radius for 
a ring thickness of 20% of the radius. Figure [D (right) 
shows the fit obtained of a compact ring, an extended 
Gaussian halo, and a 10% point source to the wavelength- 
averaged (2. 147-2. 194 /im) continuum visibilities. The fit- 
ted FWHM diameter (inner ring diameter) of the dominant 
compact model ring is 4.47 ± 0.20 mas (corresponding to 
1.12 ± 0.05 AU for a distance of 250 pc; see Tabled]). 

The measured continuum ring radius of 0.56 AU is 
~5.4 times smaller than the ^3 AU dust sublimation ra- 
dius of MWC 297 predicted by the size-luminosity re- 
lation, if there is n o shielding o f the s tellar radiation 
(jMonnier et al.l 120051: IKraus et all l2008al) . This surpris- 
ingly co mpact size was already discussed by severa l other 



authors (Eisner et al.l l2004at IMonnier et al.l I2005L [200' 



iMalbet et al.ll2007t IKraus et al.1 l2008at lAcke et al.ll200 



The correct interpretation of the compact continuum- 
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Table 1. Log of the VLTI/ AMBER/FINITO observations of MWC 297. 



Date 


Time [UT] 


AT array 


Spectral 


Wavelength DIT™ 




Calibrator 


N£ Uniform-disk 




Start End 




mode 


Rrange 






diameter of the 










H [s] 






calibrator [mas] 


2008 Apr. 06 


08:14 08:57 


E0-G0-H0 


HR-K-F d 


2.147-2.194 8 


200 


HD 175 583 


80 1.02 e 



Notes - a Detector integration time per interferogram; b Number of MWC 297 interferograms; c Number of calib rator 

interferograms; d HR mode in the K band using the fringe tracker FINITO; e Uniform-disk diameter taken from lMerand et al] 
(f2005l V 



MWC297, continuum visibility at 2.160 micron 
Two-Gaussian fit: FWHM 4.63 mas & 21 .95 mas 




MWC297, continuum visibility at 2.160 micron 
Ring diam. 4.47 mas & Gaussian FWHM 19.37 mas 



20 40 60 80 100 120 140 160 
cycles/arcsec 



20 40 60 80 100 120 140 160 
cycles/arcsec 



Fig. 2. Left: Fit of a two-component Gaussian plus a point source with a 10% flux contribution to the wavelength- 
averaged (2. 147-2. 194 /im) i-T-band continuum visibilities near Br7. Right: Fit of a compact ring plus Gaussian halo plus 
a point source with a 10% flux contribution to the iGband continuum visibilities. 
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Fig. 3. Left: Gaussian FWHM fit diameters of the Brj line-emitting region for several individual spectral channels (in 
the center and the line wings) across the Br7 emission line (red bullets). The diameters are derived from continuum- 
compensated visibilities. Therefore, they are the diameters of the Br7 line-emitting region and not of the emission region 
of both line plus continuum. The AMBER spectrum of MWC 297 is shown for comparison (green curve). Right: Gaussian 
fit to the continuum-compensated AMBER visibilities at the peak intensity of the Bi'7 emission line. 



emitted region of MWC 297 and several other Herbig Be 
stars remains disputed. Possible explanations include the 
absorption of the stellar light by an in ner gas disk allowing 
dust t o survive closer to th e star (e.g.. iMillan-Gabet et aLI 
I200H iMonnier et ail l2005t iMalbet et a l.l l2007ft. emission 
from an optically thick in ner gas disk (e.g.. IMonnier et al.l 
[2^lKraus"et~al]l2008bh . and emission from special types 
of dust grains that c an survive at higher temperatures (e.g., 
iBenistv et alll2010h . 



3.2. The characteristic size of the Brj-emitting region 

The high spectral resolution of the AMBER data allows 
us to measure the characteristic size of the Br7 line- 
emitting region for several individual spectral channels 
across the Br7 emission line (see Fig. [3] left). These diame- 
ters were determined by fitting Gaussians to the continuum- 
compensated Br7 visibilities (which are the visibilities free 
from contributions from the underlying continuum within 
the line region). The right panel of Fig. [3] shows such a 
Gaussian fit for the central wavelength of the Br7 emission 
line. 
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Table 2. Geometric model fit radii of MWC 297 in the continuum and the center of the Br7 emission line. 



Wavelength 


component CI 


diameter of CI 


component C2 


diameter of C2 


flux ratio 


range 


(dominant compact core) 




(halo) 




(C1/C2) 


Continuum" 


Gaussian 


4.63 ± 0.21 mas 


Gaussian 


£ 22 mas c 


4.6 


Continuum" 


Ring (20%) 6 


4.47 ± 0.20 mas 


Gaussian 


£ 19 mas c 


3.6 


Br7 center 


Gaussian 


12.6 ± 0.75 mas d 









Notes - a average over continuum wavelengths between 2.147 and 2.194 /im (except the line region). 6 Inner radius of a ring 
model with a ring width of 20% of the inner ring radius. c Due to the lack of data at small spatial frequencies, the size of the 
extended halo component is not well constrained by our measurements. Therefore, we can only give a rough estimate: 
22(19) ± 5 mas or larger, but smaller than the AMBER AT FOV of 250 mas. d Continuum-compensated diameter (see text). 



The continuum-compensated visibilities required for the 
size determination of the line-emitting region discussed 
above were calculated in the following way. Within the 
wavelength region of Br7, the measured visibility has 
two constituents, a pure line-emitting component and a 
continuum-emitting one, the second of which includes emis- 
sion of both the circumstellar environment and the unre- 
sol ved star. The emissio n line visibility Vb^ can be written 
as (jWeigelt et al.ll2007l) 



= v/^tot^totl 2 + \FoVo[ 2 -2F tot V tot F C V C • cos$, (1) 

where V to t (F ta t) denotes the measured total visibility (flux) 
in the Br7 line, V c (F c ) is the measured visibility (flux) in 
the continuum, and $ describes the measured wavelength- 
differential phase within the Br7 line. ^From the AMBER 
spectrum (see Fig.[T]), we find F tot = 4.5 and Fe r7 = 3.5xF c 
at the emission line center. 

We derived a large continuum-compensated radius 
(HWHM) of 6.3 ± 0.4 mas (1.6 ± 0.1 AU) for the line- 
emitting region at the peak wavelength of the emission line, 
which is ~3 times larger than the 0.56 AU HWHM radius of 
the compact continuum-emitting Gaussian co mponent. The 
smaller line-continuum size ratio obtained bv iKraus et al.l 
(2008a) can be explained by the lower spectral resolution of 
R = 1500, which leads to an averaging over different diam- 
eters corresponding to different wavelengths. Figure [3] also 
shows that the size of the line-emitting region is smaller in 
the wings of the Br7 emission line than at the line center. 
This wavelength dependence is one of several observational 
results that have to be explained by models, such as the 
disk-wind model presented in the next section. 



4. Comparison of the observations with a model of 
a magneto-centrifugally driven disk wind 

We now present a disk-wind model and compare the ob- 
served spectrum, visibilities, and phases of MWC 297 with 
the model predictions. In this disk-wind model, mass ejec- 
tion is driven by magneto-centrifugal fo rces. The theory of 
these winds was originally proposed bv lBlandford fc Pavnd 
(|1982f) and has been elabora ted by many authors , 
mainly for T Tau r i stars (see iKonigl fe Pudritzl [2000, 
iPudritz fc Baneriee" I2005L iFerreiral l2007l and references 
therein). An important property of the launching mech- 
anism in magneto-centrifugally driven disk winds is the 
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Fig. 4. Sketch of the MWC 297 disk-wind model adopted 
in this paper. 



topology of the magnetic field in the disk: magnetic field 
lines in the wind-launching region are inclined with re- 
spect to the disk plane by less than 60° (Blandford & 
Payne 1982). The radiation pressure of hot stars proba- 
bly changes the tra jectories of the gas stream s leading to 
a flatter disk wind Jde Kool fc Bcgclman 199*51 iDrew et al.l 



1998; lEverett et al.ll200lD . Figure l4l shows a sketch of the 
disk-wind model that we consider. It assumes that both an 
accretion disk and a disk wind contribute to the observed 



emission. 



4.1. Geometry and kinematics of the disk wind 

To mod el the extended disk wind, we follow an approach 
used bvlShi osman fc Vitellol ( 1993) in their studies of disk 
winds in cataclysmic variables. This approach provides a 
simple parametrization of a disk wind, which has geomet- 
rical and kinem a tical properties similar to those used by 
iKurosawa et al.l (|2006l ) for T Tauri stars (see Fig. IA.1I in 
the Appendix). Here, we use a similar description of the 
kinematical parameters of the wind where according to the 
conservation of angular momentum, the tangential velocity 
u decreases along a streamline from the Keplerian velocity 
at the base of streamline. The radial velocity v increases 
along the streamline from the initial value vq to Woo (see 
Eqs. IA.1I and IA.2I and Fig. IA.2I in Appendix A), where Vq 
and Uqo are the model parameters. 

The local mass-loss rate per unit area of the disk, m w , 
is a function of the cylindric radius cu. To describe fn w , we 
use a simple power-law 



m w (w) 



(2) 
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where 7 is a free parameter. The other model parameter is 
the total mass-loss rate 



M w = 2 



rh w (to) 2tt ujduj, 



(3) 



where uj\ and ljn are the inner and outer boundaries of the 
disk-wind ejection region, and the factor two accounts for 
the mass loss from both sides of the disk surface. 



4.2. The parameters of the accretion disk 

In classical accretion-disk models, the surface density in- 
creases toward the star as 



E(w) 



(4) 



where p = 1-1.5 (see, e.g.. lCalvet et all 120001) . At a mass 
accretion rate of M acc > 10 Moyr -1 , this disk is op- 
tically thick (due to the gas opacity) between the dust 
evaporation radius an d either the stellar magnetosphere 
(|Muzerolle et all2004j ) or the stellar surface. For MWC 297 
the v alue of M acc is fairly uncertain (see iMalbet et al.l 
|2007|) . If the surface density in the accretion disk of 
MWC 297 changes as w" 1 , then this disk will be optically 
thick up to the stellar surface. However, if the accretion 
disk extended to the innermost vicinity of the star, the disk 
luminosity at near-infrared wavelengths would exceed the 
observed one. This is the main reason why an inner disk 

gap is re quired. 

As in lAcke et all (|2008l) . we used a blackbody approxi- 
mation with a power-law distribution of the disk tempera- 
ture Trl 



Td(R) — Td(Rin){R/ Rin) a 



(5) 



where Ri n is the inner radius of the disk, and a is a parame- 
ter. The parameter a is equal to —3/4 in the two important 
special cases when a fl at disk is just heate d by the stellar 
radiation (passive disk; I Adams &: Shulll986|) or radiates du e 
to the viscous dissipation ([Lvnden- Bell fc Pringiail97l . 
However, this disk (as well as other simple power-law disk 
models) fails to rep roduce the shape of the observed con- 
tinuum visibilities (|Acke et al.ll2008l ). More sophisticated 
models can be obtained with two-component temperature- 
gradient models with slightly different temperature gradi- 
ents (see below). 



4.3. Model calculations 

For simplicity, it is assumed that the disk wind con- 
tains only hydrogen atoms. For the calculations of the 
ionization state and the number densities of the atomic 
levels, we adopted the n ume rical codes developed by 
iGrinin fc Mitskevichl (jl990l) and iTambovtseva et all (120011) 



for m oving media. These codes are based on the ISobolevI 
(1960) approximation in combination with the exact in- 
tegration of the line intensities (see Appendix A). This 
method takes into account the radiative coupling in the 
local environment of each point caused by multiple scatter- 
ing. 

The stellar rad iation of MW C 297 is described by a 
model atmosphere (|KuruczHl979l ) with an effective temper- 
ature T c ff = 25 000 K, log g = 4, and the stellar radius i?* = 



6.1 i?© as in lDrew et al.l (|l997l ). The calculations were per- 
formed for an isothermal wind. We adopted electron tem- 
perature T e values in the range 8000-10000 K, which are 
typical of the disk winds of young stars (see, e.g., Kurosawa 
et al. 2006). In our models, we adopted the fixed param- 
eters stellar mass M* = 10 M Q and a distance of 250 pc 
(iDrew et all 1 1991 . 

Figure [5] presents examples of the intensity distribu- 
tions (projected on the plane of the sky) of our best-fit 
continuum-disk and disk-wind model number 5 (see Tables 
IA.1I and l3j). The inclination angle of the model (i.e., an- 
gle between the polar axis and the viewing direction) is 
i = 20°, i.e., almost pole-on. Model images at 15 differ- 
ent wavelengths (Doppler shifts) across the Br7 line pro- 
file are shown. The images are highly asymmetric and 
wavelength-dependent. Their wavelength-dependent photo- 
center shift (with respect to the continuum image) and 
asymmetric shape change across the line, resulting in 
wavelength-dependent visibilities, differential phases, and 
closure phases. This asymmetry of the model images is the 
result of the Doppler shift of the radiation due to the gas 
motion in the emitting region. The main part of the Br7 
line emission is formed near the base of the disk wind where 
the dominant motion is Keplerian rotation. For the calcu- 
lation of the model images shown in Fig. [5j a clockwise 
motion of the disk wind was assumed. When comparing 
with the observations, we calculated the model images for 
both clockwise and anti-clockwise motion. 

Examples of theoretical Br7 line profiles of model 5 for 
different inclination angles i between the polar disk axis 
and the viewing direction (10, 20, 40, and 60°) are shown 
in Fig. [6l They suggest that a small inclination of the cir- 
cumstellar disk (20 ± 10°) is needed to explain the observed 
line profiles (Fig. 1). 

In our model calculations, we calculated more 
than 100 different models to examine the parameter 
space described in Table [31 Five examples of these 
models are listed in Table IA.1I Closest agreement is 
obtained for a two-component disk model (with two 
temperature-gradient regimes) consisting of two zones 
from i? m =0.3AU to i? s = 0.9AU and from i? s = 0.9AU 
to Rout = 3 AU separated at R s — 0.9 AU with inner and 
outer radii Ri n — 0.3 AU and R ou t = 3 AU, respectively, and 
a power-law temperature distribution in each zone given by 

Zone 1: T d = T m (R/R m ) ai ; T ln = 1800 K; a t = -0.5; 

Zone 2: T d =T m {R s /R rn )^{R/R s ) a2 ; a 2 = -0.33. 

If we choose, for example, a stronger temperature gradient, 
then the effective average disk size becomes smaller leading 
to model visibilities that are higher than those observed. 
The calculations show that a similar agreement between 
model and observations can also be obtained if the disk 
parameters vary within the ranges given by Ri n = 0.4 ± 
0.1 AU, T d (R m ) = 1700 ± 100 K, and ai =0.5 ± 0.05 up to 
R = 0.9 AU and a 2 = 0.4 ± 0.05 for R>1 AU. 

The calculations of the intensity distributions of the 
Br7-emitting region and corresponding interferometric 
model quantities (i.e., visibilities, wavelength-differential 
phases, and closure phases of the model images) show that 
they are very sensitive to the geometrical and kinematical 
parameters of the disk-wind model. Furthermore, they also 
depend on the disk inclination. We used the following ap- 
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Fig. 5. Intensity distributions of our best- fit disk- wind 
model 5 (i.e., intensity distribution of the continuum disk 
plus the disk wind; the central star is not shown; see 
Tables [3] and lA.ip at the center of the Br7 line (v = 
kms -1 ) and at 14 other velocities (the labels give the 
velocity in kms -1 ). For the calculation of the model im- 
ages in this figure, a clockwise motion of the disk wind was 
assumed. Therefore, in the blue-shifted images (left pan- 
els), mainly the disk regions on the right hand side of the 
star are bright. The radius of the inner edge of the disk 
wind ejection region (i.e., radius of the inner hole) is cji = 
17.5 i?* (^0.3 AU). The inclination angle (angle between 
the polar axis and the viewing direction) of the model is i 
= 20° (i.e., almost pole-on). The colors represent the inten- 
sity in erg ster -1 s _1 A -1 cm -2 . In these images, AMBER's 
spectral resolution of 12 000 was modeled, as described in 
Sect. \KM 




-100 100 200 

V, km/s 

Fig. 6. Normalized Br7 line profiles of the disk-wind model 
5 from Table A.l at the inclination angles i = 10° (dotted 
line), 20° (solid), 40° (dashed), and 60° (dashed-dotted; 
high spectral resolution, not degraded to the spectral reso- 
lution of AMBER) . 

Table 3. Range of parameter variations for our continuum- 
disk plus disk-wind model calculations 



Parameters 


Range 


Model 5 


Disk: 








0.25-3 AU (8.8-105 R,) 


0.3 AU (10.5 R») 


Rout 


1-5 AU (35-175 R») 


3 AU (105 R*) 


Rs 


0.85-1.25 AU (30-44 R t ) 


0.9 AU (31.5 R*) 


ai 


-0.4- -0.75 


-0.5 


OL2 


-0.34- -0.4 


-0.33 


Ti n 


1400-2000 K 


1800 K 


Disk wind: 






LJl 


0.1-3 AU (3.5-105 


0.5 AU (17.5 R») 


0J N 


0.5-5.7 AU (17.5-200 R.) 


1 AU (35 R») 


7 


-1-5 


2 


/ 


0.5-3 


0.5-3 


P 


0.3-2 


1 


6i 


10°-80° 


80° 


M w 


10" 9 -10- 6 M Q yr- 1 


10- 7 M Q yr- 1 



proach to find models that can approximately reproduce 
the observations: 

1) To select the most appropriate the disk parameters, 
we calculate the continuum intensity of the disk, Id, at 
wavelengths near Br7 and compare it with the stellar in- 
tensity /*. A good solution has to reproduce the observed 
ratio Id/ 1*, the observed if -band flux, and the observed 
visibilities. 

2) We determine the velocity and density distribution 
along each streamline corresponding to the wind geometry 
and kinematics parameters (see Sect. 2] and Appendix A). 

3) Using the velocity and density distribution, and 
choosing the temperature law along each streamline (we 
use a constant electron temperature of 8000 K), we solve 
the equations of the statistical equilibrium for the hydrogen 
atoms and compute the population of the hydrogen atomic 
levels and the ionization degree along each streamline. 

4) We then calculate the Br7 line profile and compare 
it with the observed line profile. The solution is good if 
the theoretical line profile is in good agreement with the 
intensity and shape of the observed one. 

5) For each disk-continuum and disk-wind model (see 
Table [3]) , we compute the two-dimensional intensity distri- 
bution map for several inclination angles. We calculate all 
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model intensity distributions for 15 different wavelengths 
across the Br7 line and the continuum. The spectral resolu- 
tion of AMBER of R = 12 000 is modeled in the final model 
intensity distributions (see discussion on the importance of 
modeling AMBER's spectral resolution by convolution in 
Section [O] and Fig. |A~3|) . 

6) Finally, we compute visibilities, wavelength- 
differential phases, and closure phases of the obtained 2- 
D model intensity distributions, which are needed for the 
comparison with the observations. 

4.4. Comparison of the models with the observations 

Using this method, we found that model 5 (see Tables IA.1I 
and [3]) provides the closest agreement with the AMBER 
data (see Figs. IA.41 and IA.5|) and the if-band contin- 
uum flux, as discussed below. The inner radius of the disk- 
wind ejection region of model 5 is 0.5 AU (u>x = 17.5 r*) 
and the half-opening angle of the disk wind is 80°. The 
disk-wind half-opening angle is denned as the angle be- 
tween the rotation axis and the innermost wind streamline 
(see Fig. IA.1|) . Figure [Jj shows flux, visibilities, wavelength- 
differential phases, and closure phases of model 5 for the 
same interferometric baselines as the AMBER observations 
presented in Fig. [T] 

The spectral energy distribution of model 5 predicts 
that logAF\[erg cm~ 2 s _1 ] = —7.13 in the i-T-band con- 
tinuum near Br7 for a distance of 250 pc, which is in good 
agreement with the dereddened SEP of MW C 297 from 
lAcke et al.l ([20081 ) and lAlonso-Albi et all (|2009l) . The model 
ratio Id/ 1* = 10 also agrees with the observations. 

Since the position angle of MWC 297's projected po- 
lar disk axis on the sky is unknown, the interferomet- 
ric observables of the best-fit model 5 (spectrum, visibil- 
ities, wavelength-differential phases, and closure phases) 
were calculated for all PAs in steps of 10° (see Sect. IA.3I) . 
Figure [7] compares the observations (left panel) with two 
of the models that are approximately in agreement with 
all observables (clockwise motion, i = 20°, middle: PA = 
65°, right: PA = 300°). A similar agreement can also be 
obtained for some smaller inclination angles (see Fig. [5]). 

The detailed PA dependence is presented in Sect. IA.3I 
and Figs. IA.4l and lA.5l which show the following results. For 
models with clockwise motion of the disk wind and % ~ 20° 
(i.e., almost pole-on), we obtained the closest agreement be- 
tween model and observations for the PAs of approximately 
65, 120, 245, and 300° (±5°; see Fig. EH). For the models 
with anti-clockwise motion and i ~ 20°, we obtained the 
best agreement for the PAs of approximately 15, 70, 195, 
and 250° (see Fig. lA~5|) . 

5. Discussion 

Our disk-wind model allows us to approximately repro- 
duce all interferometric observables, including the remark- 
ably strong differential and closure phases measured in our 
AMBER/VLTI observations of MWC 297 and to draw sev- 
eral conclusions: 

1) The AMBER continuum visibilities confirm previous 
resul t s dMonnier et alJl2005t iMalbet et aT1l2007t lAcke et al.l 
2008; iKraus et al.ll2008al) that the continuum-emitting re- 
gion is remarkably compact. We used geometric models 
consisting of a stellar point source, an extended halo, and 



a dominant compact Gaussian or ring to characterize its 
size. For the compact ring component, we obtained an in- 
ner ring-fit radius of ^2.2 mas (^0.56 AU). This compact 
ring radius is about 5.4 times smaller than the 3 AU dust 
sublimation radius expected for silicate grains, if there is no 
radiation-shielding material between the star and the dust 
rim. 

2) The Br7 emission line region is far more extended 
than the compact continuum-emitting region. The strong 
and wavelength-dependent closure phase reveals that the 
line-emitting region is very asymmetric along our mea- 
sured interferometric position angle and the asymmetry 
dramatically changes across the Br7 line. At the center of 
the Br7 line, we derived a Gaussian fit HWHM radius of 
~6.8 mas (1.6 AU), which is ^2.7 times larger than the 
compact 0.56 AU continuum Gaussian. The diameters of 
the Br7-line-emitting region shown in Fig. [3] are derived 
from continuum-compensated visibilities. Therefore, they 
are the diameters of the Br7-line-emitting region (not of 
the emission region of both line and continuum) . 

3) To interpret our AMBER observations, we em- 
ployed a magneto-centrifugally driven disk-wind model. 
This model consists of an accretion disk, which emits the 
observed continuum radiation, and a disk wind, which emits 
the Br7 line. The disk wind starts at a distance of ~0.5 AU 
from the star, which is further out than the typical value 
of 0.07 AU suggested for the inner ra dius of T Tauri star 
disk winds (|Safierlll993L lFerreirall2007D . The required large 
half-opening angle (~80°) of the wind in our model is 
much larger than the ty pical value (30-45°) conside red for 
T Tauri stars (see, e.g., iSafierl 119931: iFerreiral l2007h . Both 
of these features of the model are likely to be caused by 
the strong influen ce of the radiation pressure of the star 
(|Drew et al.lll998D . 

The model images of the disk wind are highly asymmet- 
ric and their photocenter shift and shape strongly change 
across the line profile, resulting in a strong wavelength de- 
pendence of the visibilities, wavelength-differential phases, 
and closure phase. This asymmetry of the model images 
is a result of the Doppler shift of the radiation caused by 
the gas motion in the emitting region. The Doppler shift 
strongly depends on the position within the disk and the 
disk inclination angle. 

The properties of the inner region of the circumstellar 
environment of MWC 297 are remarkable in many respects. 
In the NIR, it radiates at temperatures of Td < 2000 K at a 
distance of only >0.5 AU from the star. If the disk is opti- 
cally thick between the star and the sublimation zone (i.e., 
the disk is able to shield the radiation from the cental star) , 
the radius of this zone, i? S ub, will be smaller than for the 
optically thin limit. Using the model of optical ly thick and 
geometrically thin disks (|Adams fc Shul 119861) heated by 
the central star, one obtains i? sub ~ 0.8 AU for MWC 297, 
which is comparable to the inner disk radius of the best-fit 
model 5. 

The mass-loss rate of the disk wind of the best-fit model 
5 is 10~ 7 Af Q yr~ 1 (Table lA.lj) . The theory of magneto- 
centrifugal disk win ds predicts a typical ra tio M w /M acc 
on the order of 0.1 ([Konigl fc Pudrita 12000). In this case, 
the mass accretion rate in th e disk of MWC 2 97, M acc , is 
~10 -6 M Q yr -1 . According to iMuzerolle et al.l ([2004D . such 
a disk is optically thick between the star and the sublima- 
tion zone because of the gas opacity. 
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Fig. 7. Comparison of the observations (left panel; see Fig. [T]) with the corresponding model quantities of model 5 for 
two different PAs of the model image (middle and right panel; the model quantities are calculated for AMBER's spectral 
resolution of 12 000, as discussed in Sect. IA.3I and Fig. IA.3[) . The middle and right panels show the dependence of the 
interfcromctric observables (spectrum, visibilities, wavelength-differential phases, and closure phases) of our best-fit disk- 
wind model 5 (disk-wind emitting region, continuum accretion disk, plus central star; Tablesl3land lA.fi) on the wavelength 
across the Br7 line for an inclination angle of i = 20°, clock- wise motion of the disk wind, and for two different PAs of 
the projected disk polar axis on the sky: 65° (middle) and 300° (right). Several other PAs that are also approximately 
in agreement with the observations are discussed in the text and Figs. IA.4I and IA.5I The detailed dependence of the 
interferometric observables on the PA is presented in Figs. IA.4l and IA.5I 



What type of matter radiates in this compact inner 
disk inside the dust sublimation radius? We assume that 
it is a mixture of the warm, mostly neutral molecular gas 
plus refractory (e.g. graphite) grains. Simila r properties of 
the in ner compact material were found by iBenistv et al.l 
(|2010l) for the HAeBe star HD 163296. Interaction of the 
dust grains with the stellar radiation possibly plays an im- 
portant role in the acceleration of dust and gas and the 
formation of the disk wind. 

As shown above, an inner gap in the disk of MWC 297 
is needed to explain the observations. This gap region may 
be filled with material that is either f ully transparent or 
semi-transparent in the infrared (e.g., iTannirkulam et all 



2008 ). The formation of this gap could be a result of bina- 
rity (|Artvmowicz fc Lubowl[l994h . Another possibility for 
gap formation is the interaction of stellar radiation and 
wind with the inner disk. The strong radiation pressure 
and the stellar wi nd can blow away the disk atmosphere 
(jDrew et al.l fl998) . As a result, the disk can partially or 
fully dissipate in the vicinity of the star. 

4) The critical ro t ation velocity v crit — y / (2GAf*/(3i?») 
(jMaeder fc Mevnetl [2000h of MWC 297 would be 
450 kms -1 for i?» = 6 Rq and M* = 10 M Q , as adopted 
in our paper. For our determined inclination angle i = 
20 ± 10°, the maximum value of v rot sini would be 150 ± 
80 kms -1 for v ro t < Vcrit- This is in clear contradiction 
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Fig. 8. Dependence of the interferometric quantities (visibilities, wavelength-differencial phases, and closure phases) of 
our best-fit disk-wind model 5 (see Table |A~T|) on the inclination angle i: 10, 15, and 20° from left to right for the PA 300° 
of the disk (see Fig. [7]) . In contrast to Fig. [3 the Br7 line profiles are given without convolution with the instrumental 
profile of the AMBER spectrograph for clearer visibility of the theoretical line profiles. 



with the va l ue of y rn t sinz = 350 ± 50 kms 1 obtained by 
iDrew et all (fl997h for MWC 297. 

What is the cause of this discrepancy? Is it our poor 
knowledge of either mass or radius, the inclination angle 
of 20 ± 10°, or the observed v ro t sini = 350 ± 50 kms -1 ? 
The value of i is constrained not only by our modeling 
but also by t h e rece nt in terferometri c stud ies reported by 
iMalbet et all (|2007| ) and lAcke etHI (|2008T >. Nevertheless, 
further studies of i will be required to constrain i more reli- 
ably. Furthermore, it cannot be excluded that the value of 
v rot sini = 350 ±50 kms -1 is overestimated, since it is dif- 
ficult to derive the exact rotational velocity from the width 
of the observed absorption lines because of the medium 
spectral resolution employed and the quite low SNR of the 
Drew et al. spectrum caused by the relatively faint magni- 
tude of V = 12.2. Moreover, the absorption lines in MWC 
297 might be: a) contaminated by emission from outflowing 
material, as is often the case in massive young stars, or b) 
formed at a certain distance from the photosphere, where 
the broadening of the spectral line would be dominated by 
the outflow kinematics. Nevertheless, this question is far 



from being solved, and we encourage independent deter- 
minations of v rot sini using, for instance, high-resolution 
spectroscopy and state-of-the-art atmospheric modeling. 
Future multi-wavelength observations of MWC 297 using 
high spectral and spatial resolutions, combined with non- 
LTE radiative transfer modeling of the stellar atmosphere, 
will be essential to improve our knowledge of the stellar 
parameters such as the rotational velocity. 

5) According to our data, the Br7 line has a low red- 
shift in the LSR frame (but not in the heliocentric system) , 
in contrast to the model spectrum. However, it is not yet 
known whether the Br7~ emitting region is really redshifted 
with respect to the star, since the small observed redshift 
could simply be caused by the (unknown) radial velocity of 
the star (see Appendix D). 

The model of the emitting region of MWC 297 proposed 
in this paper needs to be improved. In particular, we need 
to investigate the dynamics of the magneto-centrifugal disk 
wind in the presence of the strong radiation field of the star. 
It is also necessary to investigate the effect on the thermal 
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regime of the wind of both stellar radiation and ambipolar 
diffusion. 

Additional information about the properties of the emit- 
ting region can be obtained by analyzing other hydrogen 
lines, first of all the Ha line. The observations by Drew et 
al. (1997) show that in the spectrum of MWC 297, this line 
is broader than the Bra line and has symmetric, extended 
wings up to 400-500 kms . The Br7 line is also much nar- 
rower than Ha. Our preliminary test calculations show that 
this difference can be the result of opacity effects (r (Ha) 
S> t (Br7)) in addition to Stark broadening. A more de- 
tailed discussion of the Ha line properties as well as those 
of other Balmer lines will be presented in a forthcoming 
paper. 

6. Conclusions 

We have presented VLTI/ AMBER observations of the 
Herbig Be star MWC 297 with a spectral resolution of 
12 000. This high spectral resolution has allowed us to 
study the structure and kinematics of the inner sub-AU 
accretion-ejection region in approximately ~10 different 
spectral channels (radial velocities) across the resolved B1-7 
emission line and compare the observations with disk-wind 
models. 

The if-band continuum observations confirm that 
MWC 297's continuum emitting region is remarkably com- 
pact. We measured a ring-fit radius of ~2.2 mas (—0.56 AU) 
for the dominant compact component. At all Doppler ve- 
locities across the Br7 line, the line-emitting region is both 
asymmetric (as indicated by the strong closure phase) and 
more extended than the compact continuum-emitting re- 
gion. At the center of the Br7 line, we derived a Gaussian 
continuum-compensated radius (HWHM) of 6.3 ±0.4 mas 
(1.6 ±0.1 AU) for the line-emitting region, which is —3 
times larger than the 0.58 AU HWHM of the compact 
continuum-emitting Gaussian component. 

To interpret our AMBER observations, we employed 
a magneto-centrifugally driven disk-wind model consist- 
ing of an accretion disk, which emits the observed con- 
tinuum radiation, and a disk wind, which emits the Br7 
line. We have demonstrated that all observables, i.e., the 
observed Br7 line profile, i-T-band continuum flux, visibili- 
ties, wavelength-differential phases, and remarkably strong 
closure phases, can simultaneously be reproduced by our 
disk-wind model. The modeling also shows that disk-wind 
models developed for T Tauri stars can be used to inter- 
pret Herbig Be star observations if suitable modifications 
are applied. The modification has to take into account the 
higher radiation pressure from hot stars, which can change 
the inclination angle of the streamlines of the disk wind and 
cause a flatter disk wind than in T Tauri stars. The disk 
wind of MWC 297 starts at a distance of -0.5 AU (-17.5 
stellar radii) from the star, which is a larger distance than 
the typical value of 0.07 AU predicted for T Tauri stars. 
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Appendix A: Description of the model calculations 

A.l. Disk wind 

For the description of the disk wind, we use a coordinate 
system (I, 9, w) centered at poi nt S located on the rota- 
tion axis (see the sketch in Fig. IA.1[) . In MHD models of 
the disk wind, the inclination angle of the first streamline 
with respec t to the disk plane is typ ically assumed to be 
around 60° (jBlandford fc Pavnell982|) . However, in the case 
of Herbig Be stars, the strong stellar radiation pressure pos- 
sibly bends down the streamlines and makes the disk wind 
flatter. Therefore, in our models we use a smaller inclina- 
tion angle for the innermost streamline (e.g., between 10° 
and 30°). This effect influences the characteristic size of the 
wind-launching region. 

As mentioned in Section 4, the tangential velocity of the 
wind changes along streamlines as 



(A.l) 



where uj = I sin# is the distance of the point (1, 9) from the 
rotation axis, u>i = h sin#, and uk (u>i) = (G M^/uji) 1 / 2 at 
the point (uji) at the base of streamline i. 

The radial velocity v increases along the streamlines as 



v(l) = v Q + (Voc - w ) (1 - h/lf 



(A.2) 



where vq and Voo are the initial and terminal values of the 
radial velocity and /3 is a parameter. 

We adopt = / UK^i), where UK^i) is the Keplerian 
velocity at distance uji from the disk axis and / is the 
scale factor of the asymptotic terminal velocity to the lo- 
cal Keplerian velocit y at the wind emerging point (see 
iKurosawa et al.l 120 06). The parameter / can change with 
distance w because the magnetic field in the accretion disk 
decreases with ui. An example of the velocity field of one of 
the models considered below is shown in Fig. IA.2I 

The kincmatical model described above is typ- 
ical of disk-wind models adopted for different as- 
trophysical objects with accreti on disks including T 
Tauri stars (IKurosawa et al.l [20061. cataclysmic variables 
(IShlosman fc VitelldHQQ^T and even active galactic nuclei 
dMurrav fc Chiand ll997). In these models, the gas streams 
that start in the innermost regions of the disk-wind launch- 
ing region have the highest radial velocities. These kinemat- 
ical models are characterized by non-local radiative cou- 
pling between distant parts of the moving media (s ee, e.g., 
iGrachev fc Grininl 119751: iRvbicki fc Hummer! 1 1 9781 ). which 
complicate the radiative transfer problem. In our calcula- 
tions, we neglect the influence of the non-local radiative 
coupling since its effect on the radiative excitation in the 
flattened emitting regions (such as accretion disks or disk 
winds with large opening angles) is small. 

For the calculation of the density distribution p in the 
disk wind, we use the mass-loss rate per unit solid angle, 
M w (9)/Att. One can show that M w {9) and the mass- loss 
rate per unit area of the disk, related by 



M w {9) = 4ird 2 rh w (uj) cos 



(A.3) 
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Obs. 




Fig. A.l. Sketch of the geometry of the disk- wind model adopted. 





1/ R, u/R. 

Fig. A. 2. Left: Radial velocity of the outflow along the different streamlines as a function of I; numbers 1-6 denote the 
number of the streamline. Right: The same for the rotational velocity; the dashed line shows the Keplerian velocity in 
the disk. 



where d is the distance between point S and the star. The 
total mass-loss rate is then 



M„ 



M w {9) sin 9d9. 



(A.4) 



In this case, we can write the continuity equation in its 
usual form 



Anp(l,9)v(l,9)l 2 = M w {9). 



(A.5) 



Using this equation, one can calculate the distribution of 
the number density at each point in the disk wind. 



k is the integrated line opacity in the considered spectral 
line, t(u, r) is the line optical depth at point r and for the 
frequency v in the direction to an observer, and is the 
number density of atoms in the i-th state 



t(v, r) = 



(A. 



where <f> is the profile function normalized to unity, r' is 
the vector with coordinates (x,y,z'), and z' changes from 
the current value z in Eq. IA.6I to the outer boundary value 



A. 2. Calculation of the model intensity distributions 

For the calculations of the 2-D intensity distributions 
(maps) of the emitting region, we use the coordinate sys- 
tem (x, y, z) centered on the star. The (x, y) plane coincides 
with the sky plane, the x axis is the intersection of the disk 
plane with the sky plane, and the z axis is parallel to the 
line-of-sight. In this case, the intensity of radiation at a 
frequency v within a spectral line is 

I w {v,x,y) = r m °^(r)^-i,o^)e~ T ^ r) K (r)dz,(A.6) 



where r is a vector, |r| = (x 2 + y 2 + z 2 ) 1 / 2 , v 2 (r) is the 
projection of the velocity at point r on the line-of-sight, 
and S is a source function for a transition between energy 
levels i and j 

JW^f^-lV 1 , (A.7) 



rii(r) gj 



The optical depth of the disk wind in the line frequencies 
in the direction toward the star is t*{v) = t(v, r = 0). We 
assume that there is a complete redistribution of the line 
frequencies in the reference frame of the atom and use the 
Doppler profile cj> in the calculations of the Br7 and the 
V oigt profile for the Ha line (with the same constants as 
m Kurosawa et al.ll2006t) . 

The calculations of the ionization state and the number 
densities of the atomic levels were performed in a cylinder 
with a radius r c up to 300 i?„ and a height h c = 600 i?* 
divided into grid cells in 1,9 coordinates. In each cell, we 
solved the equations of the statistical equilibrium for the 15- 
level hydrogen atom + continuum, taking into account both 
colli sion and radiative processes of excitation and ionization 
(see[G rmm fc Mitskevichl Il990l for more details). We as- 
sumed that the distribution of the atomic sub-leve ls follows 
their statistical weights. Johnson collision rates (jJohnsonl 
|1972| ) were used for all transitions except 1-2; for the lat- 
ter, we used those from lScholz et al.1 (|1990[ ). 
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Fig. A. 3. Comparison of the theoretical interferometric 
quantities of model 5 obtained with (left) and without 
(right) modeling of AMBER's spectral resolution of 12 000. 
In both cases, position and inclination angles are identical 
(PA = 300°, i = 20°). 

The main sources of continuum radiation near the Br7 
line are: i) the circumstellar disk, ii) the stellar atmosphere, 
and iii) the free-free and free-bound radiation of the disk 
wind. If the radiation of the star (/*) intersects the disk 
wind in the line of sight, it has to be corrected for the 
absorption in the wind and we observe /* exp(— t*(z/)). In 
the adopted coordinate system, we obtain for the observed 
disk intensity distribution 



I d {v, x, y) = B v (T d (R)) e 



-T w (f,x,y) 



(A.9) 



where R = (x 2 + y 2 / cos 2 i) 1 ^ 2 , T w (v,x,y) is the optical 
thickness of the disk wind at frequency v between an ob- 
server and the point in the disk with the coordinates x' = x 
and y' = y/ cosi, and i is the inclination angle (i = cor- 
responds to the pole-on orientation) . It is assumed that the 
disk is flat and optically thick; each elementary area of the 
disk radiates as a black body with the local temperature 
T d (u). 

The total intensity map of the object at frequency v is 
the sum 



I(y, %, y) = I w (v, x, y) + I d {y, x, y), 



(A.10) 



where I w (u,x,y) is the intensity of the disk wind at the 
point x, y at an emission line frequency v, and I d is the 
disk intensity distribution. 

At the central part of the map where x 2 + y 2 < 
r*, we have to add the intensity of the stellar radiation 



I if {v)exp{— (f)) to this sum. Our calculations have shown 
that the contribution of the free-bound and free-free radi- 
ation of the disk wind at the wavelength of the Br7 line is 
small compared with the disk radiation. 

To model observations with a spectral resolution lower 
than the theoretical one, the intensity maps have to be 
convolved with the instrumental profile 



Ic{v,x,y) 



/(z/, x, y)i{y' — v) dv\ 



(A.ll) 



where i(y) is the normalized instrumental profile of 
AMBER. 

In our case, the spectral resolution is quite high 
(~25 kms -1 ). The calculations of the theoretical maps 
(Fig. [5]) and the interferometric model quantities (Figs. [3 
IA.41 and IA.5[) were made using Eq. IA.11I A comparison of 
the results obtained with and without spectral convolution 
shows 



(see section IA.3|) that noticeable differences in the 
interferometric quantities exist even at this high spectral 
resolution. 

The flux at distance D from the object is an integral 
over the entire area in the x,y -plane of the emitting region 



K = J I{v, x, y)dx dy. 



(A.12) 



To summarize, the model assumptions are: 



— The disk wind is launched from the accretion disk sur- 
face in the range from distance u)\ (the starting point of 
the innermost streamline) to ujm (starting point of the 
outermost streamline). The half-opening angle (which is 
the angle between the rotation axis and the innermost 
streamline) is 9\ (see Fig. Al). 

— The wind is assumed to be isothermal with an electron 
temperature of 8000 K (for all models in Table A.l). 

— The disk wind possesses both radial and tangential (ro- 
tation) motion. 

— The number density distribution is calculated from the 
mass continuity equation for each streamline using the 
local mass-loss rate per unit area. 

— The atomic level populations were computed in the 
Sobolev approximation; a 15-level hydrogen atom has 
been considered. The 3-D radiative problem (exact in- 
tegration) has been solved in a cylinder with a radius of 
300 R* and a height 600 R* . 

A. 3. Interferometric model quantities 



Table A.l. Five examples of the computed disk- wind mod- 
els. 



Model 


M w 


9w 


Wl 


UN 


7 


/ 




W- 8 M e /yr 


deg. 


r* 


r* 






1 


3 


70 


20 


80 





3 


2 


3 


10 


30 


200 


1 


2 


3 


1.2 


80 


25 


200 


1 


2 


4 


0.6 


80 


10 


50 





* 


5 


10 


80 


17.5 


35 


2 


* 



Note: * In models 4 and 5, f is different for different streamlines 
(see Sect. A.l) 
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Fig. A. 4. Dependence of the visibilities and phases of model 5 (see Tables [3] and IA. 1|) on both the wavelength across 
the Br7 line and the PA of the model on the sky for an inclination angle of i — 20° and clockwise motion of the disk 
wind (the corresponding anti-clockwise plots are shown in Fig. IA.5|) . Top: PAs from to 80° in steps of 10° from left to 
right; Bottom: PAs = 90 to 170°. Closest agreement between observation and model is obtained for PAs of 65, 120, 245, 
a|id 300°, as discussed in Sect. IA.31 
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Fig. A. 5. Same as Fig. IA.41 but for anti-clockwise motion. The best agreement between model and observations 
obtained for the PAs of approximately 15, 70, 195, and 250° (see figure caption [A~4l) . 
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Fig. A. 6. Disk-wind model profiles of the Br7 line for different inclination angles (from top to bottom: 20, 40, and 60°). 
From left to right, the profiles correspond to models 1, 2, 3, and 4 in Table |A~T1 (model 5 is discussed in Fig. [HJ) - 



Using the theoretical intensity distribution maps of 
the emitting regions, one can calculate the interferometric 
quantities of the model: spectrum, visibilities, wavelength- 
differential phases, and closure phases for the same in- 
terferometric baseline lengths and position angles as the 
AMBER observations. Examples of interferometric quan- 
tities of disk-wind model 5 (see Table IA. II) are shown in 
Figs. \7\ E \KM \KM and IA31 



In Figure IA.3[ we show the interferometric model quan- 
tities of the Br7 line calculated for model 5 with and with- 
out convolution of the theoretical intensity maps with the 
spectral AMBER profile corresponding to spectral resolu- 
tion 12 000 (see Eg. IA.ll]) . One can see that in both cases, 
the visibilities inside the line profiles are very similar for 
all three baselines. However, the differential phases for the 
largest baseline (42 m) and the closure phases are quite 
different. 



In Figures IA.4I and IA.51 we show the detailed depen- 
dence of the visibilities and phases of model 5 (see Tables [3] 
and lA.ip on both the wavelength across the Br7 line and 
the PA of the model on the sky for both an inclination an- 
gle of i = 20° and clockwise (Fig. IA.4I) and anti-clockwise 
(Fig. IA.5[) motion of the disk wind (the PA is counted from 
north to east with north up and east to the left). In Fis. IA.41 
the closest agreement between observation and model is 
obtained for PAs of 65 and 120°. The different sign of the 
phase at the PA = 120° does not disagree with the AMBER 
observations since the sign is unknown. The plots for the 
PAs of 180 to 350° are not shown since they give the same 
results as the PAs from to 170°, except an opposite sign 
of the differential and closure phases, since a 180° rotation 
leads to a conjugate complex Fourier transform. Therefore, 
the PAs of 65 + 180 = 245° and 120 + 180 = 300° agree with 
the observations in addition to the PAs of 65 and 120°. For 
anti-clockwise motion fFig. IA.5]) . the closest agreement be- 



tween model and observations is obtained for the PAs of 
approximately 15, 70, 195, and 250°. 

A. 4. Bry line as a function of the model parameters 

Examples of the theoretical Br7 line profiles calculated for 
the models listed in Table IA.1I are presented in Fig. IA.6I It 
shows that even for large inclination angles i, the profiles 
are quite narrow. This is due to the large radius of the 
disk-wind launching region. Except for a small range of i 
values close to pole-on view (i < 20°), all orientations of the 
disk result in a double-peaked Br7 line. In observations, the 
double-peaked structure can, of course, appear single-lined 
if the line is narrow and the spectral resolution is too low. 

A single-peaked Bi - 7 emission line can also be obtained if 
a weak, narrow emission component is added to the model. 
Such emission can arise at large distances from the star, as 
discussed in the following Sect. IA.4.11 



A. 4.1. Low-velocity component of the disk wind 

A low-velocity outflow from circumstellar disks can be 
caused by the heating of the upper layers of the disk atmo- 
sphere by EUV {hv > 13.6 eV) radiation from the central 
star. As a result, the thermal velocity of atoms and ions in 
the outer part of the disk can exceed the escape velocity. 
This is the region where the disk can evaporate (Hollenbach 
et al. 1994). The inner radius of this zone is r g = GM^/vl , 
where v s is the sound speed. For ion ized gas and M* = 
10 M e , r g is -7xl0 4 /T e « 70 AU (|Gorti fc Hollenbachl 
, 2009 ). This value exceeds the radius of the AMBER FOV 
(-120 mas or -30 AU at the distance of MWC 297), mean- 
ing that the photo- evaporated disk wind cannot contribute 
to the emission observed with VLTI. However, in the ex- 
tended region of the disk between the outer radius of the 
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MHD wind ejection region (~1 AU) and r g , the EUV pho- 
tons ionize the upper layers of the disk atmosphere and 
produce a mini-H II region. This part of the disk does not 
evaporate, but can radiate in the Br7 and other lines. This 
emission is broadened by the Keplerian rotation and ther- 
mal motion of atoms. Its contribution to the total emission 
may be small, but potentially enough to fill in the small 
depression in the cental part of the double-peaked Br7 line 
profile formed by the disk wind. We simulated this profile 
(Fig. IA.7|) by changing the parameter / from / = 3 at u>\ to 
0.5 at Wat. A more detailed modeling of this narrow emis- 
sion is beyond the scope of this paper and will be presented 
in a separate one. 




-300 -200 -100 100 200 300 -300 -200 -100 100 200 300 



V, km/s V, km/s 

Fig. A. 7. Br7 line profile of model 4 with constant (a) 
and outward decreasing (b) parameter / (i = 20°; see 
Section IA.4.11 for details) . 



A. 4. 2. Screening effects 

As noted above, the optically thick part of the disk may 
screen some parts of the disk wind. This could lead to a 
red-shift of emission lines (depending on i, see Fig. IA.8|) . 
The calculations show that this effect strongly depends on 
the opening angle 9\ of the disk wind. 

At large i, the approaching part of the disk wind can 
be screened from an observer by optically thick layers of 
the flared disk. We simulated this screening effect by in- 
troducing a screen with a sharp upper boundary at some 
distance from the center. The angular size 9 of the screen 
seen from the stellar position is the parameter that we intro- 
duced. The calculations show (see Fig. IA.8|) that for some 
9 (depending on the disk-wind parameters) and i ~ 70°, 
the screening effect leads to a red-shifted emission line. It 
is likely that this mechanism is not important in our case, 
since the observations suggest that MWC 297's inclination 
angle is small (10 to 20°). 



Appendix B: Wavelength calibration 

Using the telluric spectrum provided by NOAO as a ref- 
erence (see Fig. IB.1|) . standard IRAF routines and specific 
IDL routines developed by JHG were used to find a solution 
to the pixel-to- wavelength calibration of the AMBER spec- 
trum (Fig. IB.2[) . The estimated uncertainty in the wave- 
length calibration amounts to ~2.5 kms~ 1 fFig. IB. II) . 
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Fig. A. 8. Dependence of the Br7 line profile on screening 
caused by an optically thick, dusty disk wind. The solid line 
represents the unscreened line profile of model 2 (see Table 
lA.lj) with / = 5, while the long-dashed (H d = 100 9 
= 20°) and dot-dashed (H d = 130 R*, 6 = 25°) lines show 
models with screening. For all models shown, the inclination 
angle is i = 70°. See text for further details. 

Appendix C: Beam commutation device 

To optimize the calibration accuracy, we e mployed the 
AMB ER Beam Commutation Device (BCD, iPetrov et all 
I2007D . as described in Sect. O This device is located very 
early in the optical path of the AMBER instrument and al- 
lows us to exchange two of the three telescope beams within 
a few seconds. Performing this beam commutation between 
two subsequent exposures, allows one to correct potential 
instrumental drifts. The insertion of the BCD causes a sign 
change in one of the three differential phases and the CP 
signal (|Millour et al.ll2008h . Our BCD OUT/BCD IN data 
sets precisely showed these expected sign changes. The dif- 
ferential phases and closure phases derived from both in- 
dependent data sets (BCD OUT/BCD IN) are shown in 
Fig. IC.ll Subtraction of the closure phases obtained with 
and without BCD cancels potential instrumental phase 
drifts. Fig. IC.2I presents the derived closure phase and a 
fit of a constant to the continuum closure phases. The ob- 
tained averaged continuum closure phase is —1.28 ± 0.42°. 

Appendix D: Optical spectroscopy 

MWC 297 has been rarely observed spectroscopically 
mostly because of its relative faintness (V~ 12.2 mag) and 
high reddening (B — V ~ 2.1 mag). Only two spectra that 
cover almost the entire optical range, althoug h with a low 
resolution of 0.8 A, have ever be en published (jDrew et al.l 
119971: lAndrillat fc Jaschekl Il998l) . Higher- resolution spec- 
troscopy was obtained for a few regions, mainly centered 
on strong emission or absorption lines (He I 5876 A, Nai 
D-lines, [Oi] 6300 A, and Ha, IZickgraj [200l lAcke et al.l 
2005]). A summary of the published optical spectroscopic 
data is presented in Table ID. II 

A detailed study of the object's spectrum is beyond the 
scope of this paper. Thus, we only focus on the details that 
are relevant to the modeling of our interferometric results. 
One important issue for this purpose is the geometry of 
the circumstellar environment. The strongest emission- line 
profiles exhibit a double-peaked structure, which ca n only be 
revealed at a res olving power of R > 3 000 (see IZickgraJ 
120031: lAcke et al.l l2005). This prompted [ZkkgraJ (|2003l ) to 
suggest that MWC 297 might be surrounded by a disk that 
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Fig. B.l. Wavelength calibration of the AMBER data. 
Top: Raw AMBER spectrum of MWC 297 (solid blue 
line) compared to a telluric spectrum from the Kitt Peak 
Observatory. The green line shows the original Kitt Peak 
spectrum with a spectral resolution of 60 000, while the 
red line displays the telluric spectrum that has been con- 
volved to the same spectral resolution as the AMBER data 
(R = 12 000). For clarity, the red curve has been shifted 
down by 0.4. As the figure reveals, after the wavelength cal- 
ibration, the AMBER raw spectrum of MWC 297 matches 
several telluric features. Bottom: Same as top panel, but 
just for the wavelength region around the Br 7 line. 



is viewed at an intermediate inclination with respect to the 
line-of-sight. 

Another crucial parameter is the systemic (stellar) ra- 
dial velocity. Since this had never been carefully discussed, 
we critically examined possible methods to derive this 
quantity and performed an independent study of the pub- 
lished spectra. Absorption lines of the photospheric origi n 
are only seen in the spectra obtained bv lDrew et al.l (|l997h . 
They include a few He I lines blueward of 5800 A and those 
of C in and N in in the 4500-4700 A range. All of them are 
broad, which was interpreted to be caused by a high stellar 
rotation rate (jDrew et al.lll997l) . Additionally, the residual 
intensities of the carbon and nitrogen lines are within ~5 % 
of the continuum, hampering high accuracy measurements 
even at the SNR of ~ 120 achieved in the spectrum. The av- 
erage LSR radial velocity, which we derived from five He I 
absorption lines at 4009, 4026, 4143, 4471, and 4713 A, is 
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Fig. B.2. Top: Calibration between the pixel and 
wavelength values of the VLTI/ AMBER observations of 
MWC 297. A least-square polynomial fit of the data is 
also shown (red line). Bottom: Residuals between the mea- 
sured position of the telluric line spectrum obtained from 
the wavelength calibration of the AMBER data and the 
expected position from the reference telluric spectrum pro- 
vided by NO AO. The root- mean square value is of the order 
of 2.5 kms -1 . 



— 7 ± 7 km s -1 . However, He I lines in the spectra of B-type 
emission-line stars are known to be strongly affected by the 
circumstellar material in the dense regions close to the star. 
Even pure absorption lines can exhibit variations that re- 
sult in the observed shift s from the photospheric positions 
fsee llsraelian et aT1ll996[ ). Therefore, these lines cannot be 
considered as a reliable source of the systemic velocity. The 
most clearly recognized non-helium photospheric line is the 
C in 4650.16 A line, which however, forms a blend with two 
Nm lines at 4640.64 and 4641.90 A. This line gives a more 
negative LSR radial velocity of —29 kms -1 . The accuracy 
of this value is at least of the order of 15 kms -1 . The line 
can also be affected by the circumstellar material, which 
also ensures that this measurement is unreliable. 

An alternative method for deriving the systemic ve- 
locity can be used if photospheric lines are not seen in 
the spectrum or there are no data of sufficient quality for 
these lines. It employs centroid radial velocities of emission 
lines, such as those of Fen. It has been successfully tested 
on co mplex objects, such as LBVs (e.g.. iHumphrevs et al.l 
1989) and stars that exhib it the B[e] phenomenon (e.g., 
iMiroshnichenko et alJl2002T ). The resulting systemic veloc- 
ity is used to determine a kinematical distance and lumi- 
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velocity (LSR) [km/s] velocity (LSR) [km/s] 

-100 100 -100 100 
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Fig. C.l. Differential phases and closure phases derived 
from the data obtained without (left, BCD OUT) and with 
(right, BCD IN) the beam commutation device. The shown 
signals are not yet corrected using the calibrator data. After 
we corrected the sign changes introduced by the BCD de- 
vice, the two data sets yield the shown signals, which are 
nearly identical and demonstrate the high accuracy of the 
measured phase signals. 



Closure Phase: (BCD-OUT - BCD-IN)/2 
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Fig. C.2. BCD-calibrated closure phase signal obtained by 
subtracting the calibrated closure phases obtained from the 
BCD IN and BCD OUT data. The blue line is the fit of a 
constant to the continuum (2.147-2.194/im) closure phase 
outside the B1-7 line. 

nosity of the object that are in general good agreement 
with those found by other methods (e.g., spectroscopic par- 
allax, interstellar extinction versus distance relationship). 
The method seems to work for any type of centrally sym- 
metric circumstellar envelopes. Moreover, it has been shown 
to trace the central sta r motion even in binary Be stars (e.g., 
Harmane c et aLlfeOOOh . whose disks may con tain asymmet- 
ric density_r^turbations fsee lOkazakilll99ll ). 

I Drew et all |l997) comment that the centroid LSR ra- 
dial velocities of the Ha, H/3, and the Hei 5876 A lines are 
located between +2 and +20 kms -1 . Other measurements, 
including those we made in the optical spectra and which 
we have at our disposal, are listed in Table ID.ll The av- 
erag e value is +8 ± 1 kms -1 , similar to +10 ± 2 kms -1 
from lAcke et al.l (|2008T ) . We note that the formal statistical 
uncertainty of these velocities may be larger, because the 
number of lines used is relatively small. In addition, a sys- 
tematic error due to, for instance, a slight asymm etry in the 
line profiles (see those shown by [Zickgraf 2003|) cannot be 
excluded. Therefore, we assume that the real uncertainty 
in the systemic velocity estimated with this method might 
be as large as ^5 kms -1 . 



Table D.l. Summary of published optical spectroscopic 
observations of MWC 297 



Date AA, R V Ref. 

A A/AA kms -1 

1987 Sept. 9 6535-6600 23 000 1 

1987 Sept. 13 5868-5903 45 000 +6 ± l c 1 

6290-6335 1 

1992 Sept. 9 6290-6759 8 000 a 2 

1994 June 20-22 3876-5198 8 000 a +10±8 d 3 

5840-8809 3 

1997 July 21, 23 4119-5734 5000 b 4 

1997 July 20, 24, 25 6297-8911 4 

2002 Apr. 1 6279-6321 125 000 +9 ± l c 5 

2002 June 25 5635-6637 30 000 +10±2 f 5 



Column 1: observing date; Column 2: observed spectral 
range; Column 3: average resolving power; Column 4: av- 
erage LSR radial velocity (see explanations below and in 
the text) ; Column 5: reference to the reporting paper. 

a - two-pixel resolution 0.8 A 

b - 1.3 A resolution at A = 6000 A 

c - from the [O 1] 6300 A line 

d - from emission lines of Fell 4923 and 5018 A [Fell] 5159 
A, Sill 5041 and 5056 A, and Sill 5958 and 5979 A 
- from the [O 1] 6300 and 6363 A lines 

f - result publis hed bvlAcke et al.1 (120081) 

R eferences: 1 - IZickgraJ ll2003h: 2 - IVieira et all (120031 ): 3 
- iDrew et ail (119971 ); 4 - lAndrillat fc Jaschekl (jl998l ); 5 - 
lAcke et all l|2005h 



To summarize the above discussion, we conclude that 
the positions of the emission lines suggest a systemic radial 
velocity of MWC 297 of +8 ±5 kms" 1 . In contrast, the Hei 
absorption lines give —7 kms -1 with an even larger uncer- 
tainty, because of their possible displacement from the pho- 
tospheric position due to processing by the circumstellar 
material. The value derived from the emission lines seems 
to be more reliable because it is based on data from several 
spectra obtained over a period of 15 years. Nevertheless, 
it still needs to be confirmed using high-resolution data of 
more lines. 
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